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Abstract 



Qh| We pursue the study of the high energy behaviour of the gluon propagator on 

O ; the lattice in the Landau gauge in the flavorless case (n/ = 0). It was shown in 

O , a preceding paper that the gluon propagator did not reach three-loop asymptotic 

scaling at an energy scale as high as 5 GeV. Our present high statistics analysis 

includes also a simulation at /? = 6.8 (a ~ 0.03 fm), which allows to reach /x ~ 

^ . 10 GeV. Special care has been devoted to the finite lattice-spacing artifacts as well 

H I as to the finite volume effects, the latter being acute at /3 = 6.8 where the volume 

is bounded by technical limits. Our main conclusion is a strong evidence that 

the gluon propagator has reached three-loop asymptotic scaling, at fi ranging from 

5.6 GeV to 9.5 GeV. We buttress up this conclusion on several demanding criteria 

of asymptoticity, including scheme independence. Our fit in the 5.6 GeV to 9.5 GeV 

window yields A = 319 it 14 1^2o MeV, in good agreement with our previous 

result, A = 295 it 20 MeV, obtained from the three gluon vertex, but it is 

significantly above the Schrodinger functional method estimate : 238 it 19 MeV. 

The latter difference is not understood. Confirming our previous paper, we show 

that a fourth loop is necessary to fit the whole (2.8 -i- 9.5) GeV energy window. 

LPT Orsay/99-72 



1 Unite Mixte de Recherche - UMR 8627 

^ Unite Mixte de Recherche C7644 du Centre National de la Recherche Scientifique 

e-mail : Philippe . BoucaudOth . u-psud . f r , roiesnelficpht . polytechnique . f r 



In previous works, we tackled the non-perturbative calculation of the QCD running 
coupling constant in two different ways: (i) by using the three gluon coupling [|I|, 0] , and (ii) 
by matching the behaviour of the lattice regularized gluon propagator to the one predicted 
by perturbation theory [Q . The latter method was expected to take benefit of the very 
good statistical accuracy of the propagators and thus to yield a rather precise estimate 
of the strong coupling constant in its ultraviolet (UV) regime, i.e. of Aqcd- Unluckily 
we could not fulfill this program for the unexpected reason that the gluon propagator 
has not yet reached the asymptotic scaling at scales of 2.5 -7- 5.0 GeV. This conclusion was 
supported by several different tests. In particular, the remainder of a strong scheme 
dependence when using one-, two- and three-loop formulae indicated the compelling need 
of higher orders in the perturbative expansion. Still we observed that the inclusion of 
third-loop corrections improved the asymptoticity (albeit not enough) over the two-loop 
results. We were naturally tempted to extend the analysis to higher energy scales where 
any perturbative expansion, with a fixed number of terms, should progressively improve. 
This is the basic motivation of the present paper. 

Since we want to reach ever larger momenta on the lattice, we have to assure that the 
dominant lattice artifacts are under control. We must also ensure that the energy window, 
in which we could test scaling of the gluon propagator, is large enough. The reason is 
that at large scales, the coupling constant has a very mild logarithmic dependence. If 
the window is too narrow, the higher order terms might mimic the lower order ones, 
that we consider, and thus introduce a bias in Aqcd- A very wide energy window has 
been explored in ref. [J18l by using the Schrodinger functional technique. When using the 



methods based on Green functions |I|, 0, Q, it is more difficult to vary the energy scales 
on several orders of magnitude. Moreover, one deals with more scales: the lattice spacing 
a, the linear lattice extension L, and the momenta of the gluons p"^. On the other hand, 
compared to the Schrodinger functional method, we believe that the Green functions have 
a simpler physical meaning and, being conceptually very different, represent a necessary 
test. 

The requirement of larger momentum scales implies smaller lattice spacings if we are to 
keep the UV 0{a?'p^) artifacts under control [ap -C 1). Equivalently, we need to perform 
the simulations at larger /3, which (for reasonable computing time) also means smaller 
volumes and potentially dangerous infrared (IR) finite volume artifacts. To prevent these 
problems, we need to ensure Lp ^ 1. 

The question is whether we can find an ensemble of lattice results fulfilling all 
requirements, i.e. that the lattice artifacts are small enough (0(L~^) -C p -C 0{a^^)), 
and that the energy window is large. This is a very demanding requirement because a 
small change in the coupling constant can induce a large uncertainty in Aqcd- To avoid 
the statistical uncertainties in that respect, we work with 1.000 independent gauge field 
configurations at every stage of this research R. To keep the energy window large enough. 



^By a lattice result, we mean a value of the bare propagator for one value of p^ obtained on a lattice 
with a given /?, and in a specific volume. 

^The only exception is the simulation at (6.0, 32"'), where we have 100 configurations. 



we need to fit simultaneously the lattice data obtained at different lattice spacings. In 
other words, we must ensure that the small momentum lattice data with large /3, and 
large momentum data with small (but reasonable) /?, are compatible. In order to achieve 
that, and to reduce systematic uncertainties to the level of statistical ones, one evidently 
needs to control both IR and UV lattice artifacts for all the data. Therefore, a major 
part of this paper is devoted to these issues: the reduction of the UV and IR systematic 
uncertainties. The discretization errors are monitored by working at three values of the 
bare lattice coupling: 

(3 = { 6.0, 6.2, 6.8 } . (1) 

A lattice spacing a~^ = 2.72(11) GeV has recently been measured at /5 = 6.2 with a 
non-perturbatively improved action |T^. For a direct comparison with |jl|, 0, |^, in this 



paper we will keep a ^{(3 = 6.2) = 2.75 GeV which is well within the error bars. Other 
lattices are calibrated relatively to this one, by using the lattice measurement of the string 



tension |15|. We take: 

a-i = { 1.97 GeV, 2.75 GeV, 6.10 GeV } , (2) 

at /3 = 6.0, 6.2 and 6.8 respectively. Thus, our lattice spacings vary from 0.03 fm to 0.10 fm. 
The study at /? = 6.8 (a ~ 0.03 fm) allows us to reach momenta up to ~ 10 GeV. The main 
study of the finite volume effects is performed at /? = 6.0, by repeating the calculation 
with the following lattice volumes: 

V = { l2^ l6^ 24^ 32^ } . (3) 



From this study we deduced an efficient parametrization of the finite volume effects 
which allows us to extrapolate our high momentum data to the l^ — >■ cxo limit. An 
important cross check is provided by two volumes at /3 = 6.8: 

V = { 16^ 24^ } , (4) 

while for j3 = 6.2, we work with the volume V = 24^ only. 

An improved version of the method used in ref. to cure lattice hypercubic effects, 
has been applied. The non-hyper cubic finite spacing effects have been dealt with by 
comparing different values of (3. 

Curing the above-mentioned artifacts, we could keep most of our lattice points. We 
discarded those which exhibit large corrections. 

In this paper, we will not deal with the small momentum behaviour of the gluon prop- 
agator. We postpone it to our forthcoming publications. This part has so far attracted a 
lot of attention in the literature [0-|jTO[ (for a review with a rather complete list of refer- 
ences, see ref. [0), while, to our knowledge, the large momentum gluon propagator has 
not been studied in detail. Preliminary results of this paper were presented in ref. |I2|. 



The remainder of this paper is organized as follows: In section 1 we outline the gen- 
eralities of the method we use (previously described in ref. H), and introduce the main 



perturbation theory tools. In section 2, the lattice artifacts are discussed. Those related 
to hypercubic geometry are eliminated by improving the method presented in ref. 0. A 
procedure allowing to treat empirically the finite- volume effects is described. In section 3 
we perform a three-loop fit to the gluon propagator in the MOM scheme and apply the 
test of scheme independence, by considering a (5.6 -v- 9.5) GeV window in which only the 
data at /5 = 6.8 are used. Section 4 is devoted to the study of the whole window, ranging 
from 2.8 up to 9.5 GeV by combining all lattice data. We discuss our results in section 5 
and conclude in section 6. 

1 General description of the method 

The Euclidean two point Green function in momentum space writes in the Landau gauge: 

G^^lTip, -P) = 5„,., (5,,,, - ^) G^'\p') , (5) 

where ai, 02 are the color indices ranging from 1 to 8. The bare gluon propagator in the 
Landau gauge (see for instance ref. [§), is such that 



A JT„ .,2 A^^ /Vln 1,2 ' ^ ' 



A-*oo dlnfi'^ A^oo dln/i^ 

is independent of any regularization scheme. Zs{fi,A) is the gluon renormalization con- 
stant in the MOM (or MOM) scheme at the point p"^ = fJ''^, and A is a generic notation 
for the UV cut-off (a"^ or {d - A)~^). 

It is well known that the /x and the A dependences of Z3(/x,A) factorize when one 
drops all the terms vanishing as A ^ 00 (see ref. ||13|), and we can write: 

Zsifi,A) = Z^{^)Z',{A) + Oil/A) ■ (7) 

where the evolution of both, Z^{fi) and Z|(A), is described by Callan-Symanzik equations 

^ - r^i^)) z^{^^) = , 



(iln/i^ 
d 



T\A) Zl{A) = . 



dlnA'^ 
From the QCD perturbation theory we know that 

dlnZ^ifi) r^ f-fo 71 „2 , 72 „3 , ^/ 4^ 



dlnn^ V47r (47r)2 (47r)3 J 

where it is understood that the coupling constant in a given scheme is a function of fi 
such that 

^'^ f^r \ Po 2 A3 P2 i . ^i 5n /,^n 

— p[a) = —-r-a — . a — a + 0[a ) (10) 
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with 

/3o = ll, A = 51, 70 = y, (11) 

in the flavorless case {nj = 0), while 71,72 and P2 are scheme dependent. To be specific, 
in the fiavorless MOM scheme^: 

29 

/32^4824., 71 = ^, 72^960. (12) 

o 

Lattice calculations provide us with the bare propagator but in a finite volume which, 
besides the UV cut-off (A ~ 1/a), introduces an additional length dimension, L (the 
physical volume being L'^). As we shall see, finite volume effects, as well as hypercubic 
artifacts, should be eliminated first in order to have access to the renormalization constant 
in eq. (§), Z3(/i, 1/a). 

Equations formally analogous to (§) and (|l^) can be obtained from the second line 
of eq. (|), with the substitutions of 1/a for A and of the lattice bare coupling constant, 
a^ = 3/(27r/3), for the renormalized one. Unhappily the anomalous dimension coefficients, 
7J, 72, have not been determined to our knowledge, presumably due to the difficulty of 
the task. Any perturbative calculation of ^3(1/0) appears thus to be limited to one loop, 

dlnZlJl/a) _7o^ 1 
da'il/a) /?o«'(l/a) 

for which it can easily be proven that 7q = —70- 

Our general ''strategy", as explained in ref. [Q, will be to integrate simultaneously 
eqs. (||,0), up to three loops in a given scheme. The solutions depend on the initial 
values Z|^(/io) and a(/io)- They are related to the lattice results, Z^i^fi, 1/a), through 
eq. (^. When we restrict our computation to the lattice data at /9 = 6.8, Z|(l/a) is 
just an overall irrelevant multiplicative constant. With our three values of f3, two ratios 
of three {Zl{Q.8)/Zl{6.2), Zl{6.8)/Zl{6.0), Zl{Q.2)/Zl{6.0)) are necessary for appropriate 
matching of all our lattice data. These ratios will be fitted and compared to one-loop 
predictions in section 4. 

Finally, the knowledge of a(/Uo) in a given scheme allows the determination of Aqcd 
in this scheme, and hence of Aq^j-,. 

2 Lattice artifacts 

We refer to ref. [^ for technical details concerning the lattice setup in our simulations, 
the calculation of the Green functions, their Fourier transform, the checks of the 5ai,a2 
color dependence of the propagators, and the set of momenta considered for the different 
lattices studied. Since the release of ref. |]I[, we increased the statistical quality of our 
data, and further explored various lattice volumes and various values of /3. As mentioned 

•^Details on the computation of the parameters /32, 71, 72 in this scheme can be found in refs. ^ Q. 



in the introduction, of special interest for this study are the results of our simulation 
performed at (3 = 6.8, at two volumes: 16^ and 24*^. The high statistical accuracy of our 
data made a detailed study of systematic uncertainties possible and mandatory. 



2.1 Hypercubic artifacts and other (a p ) effects 

We start with the discretization errors. In a finite hypercubic volume the momenta are 

the discrete sets of vectors 

27r 



Pm 



n 



V- ' 



(14) 



where the components of n^ are integers and L is the lattice size. The propagators 
have been averaged as usual, over the hypercubic isometry group if 4. The momenta 
corresponding to different orbits of if4, but belonging to the same orbit of the continuum 
isometry group 50(4) {e.g. n^ = (2,0,0,0), and n^ = (1,1,1,1)), have been analyzed 
according to an improved version of the method proposed in ref. [Q. 

Let us briefly recall the elements of that method. The main idea is based on the fact 
that, on the lattice, an invariant scalar form factor, like G^^\p'^), is indeed a function 
of 4 invariants p^"' = ^„p!J, n = 2,4,6,8. We will neglect the invariants with degree 
higher than 4 since, in any case, they vanish at least as ~ a"^. Thus, we parametrize 
and expand the lattice two-point scalar form factor as a function of the two remaining 
invariants which, on dimensional grounds, appear as p"^ and a^p^: 



G^lliP^a^P 



'-^'hL,a) 



G'g(p^O;L,a) + 



dG 



Mat ' 



a(a2pw; 



aV''^ 



a2p[41=0 



(15) 



dG/d{a^p 



[4^ 



2„[41 



I a^p 



^Q is a symbolic notation for the lattice hypercubic slope. This equation 



summarizes our method to reduce hypercubic artifacts and contains all the assumptions 
on which the method relies. We make the hypothesis that the lattice propagator for the 
discrete momenta belonging to different if4-orbits takes values according to a certain con- 
tinuous functional behaviour. When several orbits exist for one p"^, to the extent that a 
linear approximation is licit, the hypercubic slope can be extracted and the extrapolation 
to a^p^^' = 0, by using eq. (15), can be made. We have checked that the linear approxima- 
tion ( [T5| ) is indeed good enough. This is what has been done in ref. 0]. Now we elaborate 
further on this method in order to improve its accuracy and extend its applicability. 

A simple dimensional argument leads to: dG/{a'^dp^^^)\^2 m-n "^ '^Iv^-, as L ^ 00, 
suggesting the fitting form &/p^, where 6 is a constant |. At this point, from the study of 
the lattice hypercubic slope for different values of /?, we extracted the values of term 
hjp'' and examined the behaviour of the remainder. We find that the logarithm of 



-6 + p4aG'/(a2apW)L: 



2p[4]=0 



is well described by a linear function of Lp. This leads 



be ^ 



*For example, the hypercubic correction for the free propagator, 1/p^ = Xjij? — j^a^pW) + . . ., would 



us to: 



dG 



9(a2pW> 



a2pW=o 



p4 



f 1 + c exp(— d Lp) 



(16) 



By fitting our data to the above form and by keeping y^ /d.o.f. ~ 1 (see Figjl|a), we 
obtained the following set of parameters: 



/? 


h c d 


6.0 


0.117(9) 1.1(1.5) 0.15(8) 


6.2 


0.109(2) 1.9(1.2) 0.18(3) 


6.8 


0.100(5) 0.6(8) 0.13(5) 



Table 1: The values of the parameters h, c, d, as obtained by fitting our data to eq.(|T 



The values for c and d in Tab. |l] are presented to show the order of magnitude. Their 
errors, estimated by using the jackknife method, are misleading since the parameters are 
strongly correlated. A refined statistical study is not necessary for our purpose, since we 
follow the jackknife analysis cluster by cluster to the end. 

On the other hand it is rewarding that the errors on b in (|1|) are small, which is essential 
for an accurate infinite volume limit. It is also encouraging that b varies only slowly with 
a, which justifies our neglecting higher order terms, 0{a'^ pf^') etc., and it clearly confirms 
that we do control the lattice hypercubic artifacts. We can now take advantage of the 
good fit obtained with eq. (|l6|) and compute the lattice hypercubic slope in cases where 
only one orbit exists (for instance, when n^ = 5, we only have n^ = (2, 1, 0, 0) and its H4 
orbit). 

To summarize, by using eq. (p^Bf), we extrapolate G|J(p^, a^pW; L, a) to what we will call 
the hypercubic-free propagator, GU(p^, 0; L, a). In ref. 0, we discussed the improvement 
brought in by our previous (a^p'^' = 0)-extrapolation approach with respect to other 
methods to reduce hypercubic artifacts. The use of eq. (|TB|) allows even better accuracy 
on slopes: not only does it allow an extrapolation to a^p^^' = when only one orbit exists 
but it also helps to reduce the uncertainty by taking benefit of the neighboring values of 
n'^ when the error is locally too large. The outcome of this improvement is depicted in 
Fig. Pj, where the resulting curve joining the points obtained by using eq. ([TBI) is much 
smoother than the one obtained by separate extrapolation for each n^. 

It is important to add that not all O(a^) artifacts are eliminated: for instance, the 
lattice artifacts oc a^ p^, which do not break 5*0(4) invariance, are still present. One way 

^Note however, that even in the latter case, the resuhing curve is by far smoother than the one 
obtained by simply averaging the orbits or by selecting the "democratic" points, as advocated in ref. ||]. 
This improvement was already illustrated in ref. |p[ . 
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Figure 1: Plot (a) shows p - 



2 dG 



a2pW=0 



, as a function of the scale p evaluated on 24'^ 



lattices at (3 = 6.0, 6.2, 6.8. Plot (b) compares the "hypercubic-free propagator" defined 
from eqs. (]IJ,|I^j (black circles) and the propagator computed by direct extrapolation to 
a2p[4] _ g (^^]2ite circles), plotted as a function of the momentum. 



to deal with this problem is to compare the data at different values of (3. One can also 
study the stability of the results when the maximum value of p^ used in the fits is varied. 
We have checked this stability in the fits presented below. 

2.2 Finite volume effects 

After removing hypercubic artifacts, we are left with the hypercubic-free propagator. The 
dependence on the length scale, L, as we previously mentioned, is apparent from Fig. ^ 
The elimination of this additional length scale should be done in order to compute the 
renormalization constant for the Landau gauge gluon propagator, 

Z3(/i, 1/a) = hm ip'G[%^ 0; L, a))^.^^, . (17) 

In doing so, we will not attempt a theoretical understanding of the expected finite volume 
dependence of the bare propagator. We will be content if we obtain a reliable empirical 
parametrization for the dependence on the lattice volume of GU(p^,0;L,a) which will 
allow us to take the required limit (|T^. For dimensional reasons, we will take it as a 
function of Lp and a/L. We note that the difference amongst the data at fixed (3 and 
various volumes gets bigger as we move towards lower p^. This is illustrated in Fig. |^ where 
we plot p^GU(p^, 0; L, a). We tried several ansdtze for this parametrization function and 
finally found that 

G[ll{p^O;L,a) = ^^(p^, 0; oo, a) (l+v,(^^y + V2exp{-V3Lp) 
with 

vi = 450(40), V2 = 0.44(16), V3 = 0.177(25), (19) 

gives the best fit to the behaviour on L of hypercubic-free propagators evaluated on 
12^, 16^ 24^ 32^ lattices in the energy window f\ (1.5 ^ 3.0) GeV at /5 = 6.0. This 
parametrization is not efficient at lower energies, as it can be seen in Fig. Qb. 

Once the parametrization function from propagators evaluated at /? = 6.0 is estab- 
lished, it can be applied to our results at 16^ and 24^ lattices at /3 = 6.8. The agreement 
after extrapolation shown by the curves resulting from the extrapolation in Fig. ^ is a 
crucial test for the validity of such a parametrization for the finite volume effects. 

Thus, eqs. (0-[l^) lead to a non-perturbative evaluation of the renormalization con- 
stant, Z^lfi, 1/a), which is precisely the quantity we want to compare to the predictions 
of perturbative QCD. 

^In this procedure we assume that the volume 32'* is aheady infinite, i.e. we take 
G[ll{p'^,0;32,a{6.0)) ~ G[ll{p'^,0;oo,a{6.Q)). In the fit to the form (|l|), the energy window is cho- 
sen such that the total x^ /d.o.f. ^ 2. 
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Figure 2: Plot (a) contains the hypercubic-free propagators evaluated on 12'* (black 
squares), 16^ (black circles), 24"^ (white circles), 32"^ fwiite squares) lattices at [3 = 6.0. 
The same data extrapolated to L ^ oo according to the parametrization given by 



eqs. ( jlq , !!^ ) are plotted in (b) . 



3 Fitting in the MOM scheme at high energies 

In this section we perform the matching of the perturbative predictions for MOM renor- 
mahzation constant to our lattice result at high energies. We follow the method outlined 
in section |l| and refer the reader to ref. [^ for more details. 

We consider the coupled differential equations (§,0) in the MOM scheme, where the 



coefficients 71, 72, /32 are those given in eq. (jT^). We fit our data at /? = 6.8 for Z^^j,, 1/a), 
with a solution of these coupled equations in the energy window (5.6 -^ 9.5) GeV. The 
result of this fit is 



^3(9.5 GeV) = 1.3107(9), 



a 



MOM 



(9.5 GeV) = 0.190(3), x^/d.o.f. = 0.29 . (20) 



The x^/d-o.f is significantly smaller than 1 which may be a sign of some correlation 
between the points at different values of the energy /i. 

As explained in |^, A'^^ can be estimated from the above quantities (^), by using 
the perturbative expressions to two- and three-loop accuracy. We obtain: 



A MS 
^^(31oop) 

aMS 
-'^(21oop) 



0.346 A^?Z) 



0.346 A 



MOM 
(21oop) 



319 ± 14 MeV, 
375 MeV, 



(21) 



10 



(a) 



1.6 



^ 1.5 

N 



1.3 



^ = 6.8 

« , o 16' 

# '^ .24' 



• 



§ 



•. 5 






o o 




4.5 5.5 6.5 7.5 S.5 9.5 4.5 5.5 6.5 7.5 8.5 9.5 

\i [GeV] ]i. [GeV] 

Figure 3: This figure is analogous to Fig. ^ for (3 = 6.8. Black (white) circles correspond 
to 24^ (16^; lattices. 



where the error is only statistical at this stage. The existence of a good fit (see Fig. 1) 
by itself is not a sufficient proof of asymptoticity: next-to-three-loop corrections could be 
mimicked by a simple rescaling of Aqcd in the considered energy range f^. This is why we 
developed a consistent method to test asymptoticity by exploring the scheme dependence 
within the domain of so-called ^^good schemes": one investigates the dispersion of the 
result for A'^^ when we vary the schemes , by varying 71 and 72, in all the possible ways 
such that the successive terms in the perturbative series (p|,p!OD are at most as large as the 
preceding ones []. In ref. 0, we found that this dispersion is of ~ 35 MeV, for A^^ fixed 
from the fit of the gluon propagator at ~ 4 GeV to the three-loop perturbative expression. 
In the present study, when fixing A'^^ at around 9.5 GeV, this dispersion (in the same 
domain of schemes) is of ~ 10 MeV only. 

Note that the difference in eq. (pTD between two and three loops, although smaller 
than at 4 GeV, is still sizable indicating the necessity to include the third loop term. 

4 Global description and asymptotic pattern 

Our analysis at high energy in the previous section seems to establish that a signal of 
three-loop perturbative scaling is found in the energy window (5.6^9.5) GeV. By including 
the data obtained at /3 = 6.0 and (3 = 6.2, one can make the energy window larger: 2.8 
GeV< /i < 9.5 GeV (the choice for the lower limit will be discussed below). However, as 



^This is the generalization of the effective charge approach proposed in ref. |T3| 
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|x [GeV] 

Figure 4: The plot shows (full line) the best fit to Z^lfi, 1/a) with the three-loop formula 
together with the lattice results at j3 = 6.8 after extrapolation to L — > cx) from both 16*^ 
and 24^ lattice volumes, for 5.6 GeV < /x < 9.5GeV. The fit is continued outside the 
energy window as a dashed line. 



already mentioned in the end of Sec. |T], any global fit involves two additional parameters, 
the ratios Z|(6. 0)7^3(6. 8) and Z|(6.2)/Z|(6.8). Moreover, we do not expect a three-loop 
perturbative behaviour to work in the whole energy window. In ref. 0, we showed that 
the propagator was not asymptotic to three-loops at ~ 4 GeV. The difference between 
^(3k)op) ^^ ^'i- (HI) ^^^ ~ 3^'-' MeV as found in ref. P], confirms that statement. Thus, 
at least the fourth loop correction is necessary for the global fit. Unfortunately, such a 

perturbative result is not available in the MOM scheme. 

With these fives free parameters, a global fit turns out to be unstable. A global study 
of our lattice data would nevertheless enable a direct test of consistency for the the whole 
information we extract from gluon propagator. For that reason we adopt the following 
strategy. First, we take the value of AM|^ -^ given in eq. (0), to be the asymptotic one, 
i.e. we assume the gluon propagator to reach asymptotia at three-loop for the energy 
window studied in the previous section. Then we fix the fourth loop correction to the 
three-loop perturbative expression, by fitting the data obtained in our simulations at 
f3 = 6.0 and f3 = 6.2, corresponding to the energy range (2.8 -h 4.3) GeV. Once, the fourth 
loop correction is known, we will verify the asymptoticity of the gluon propagator in the 
entire energy window (2.8 -7- 9.5) GeV. 

The four-loop information about the gluon propagator in MOM-scheme is encoded in 
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the coefficients 7!^'-'^ and /3|^°^. These two coefficients are not independent but related 
through the expression |1^, |^ 

(47r)4^ (47r)4/5o (47r)4 ^ (47r) ^o ' ^ ^ 

which is vahd for any renormahzation scheme in which 71 = ^yj^o^^^ 72 = ^^"-"^ and 

P2 = ^2^'^^, hsted in eq. (|T2|). Thus, there is only one free parameter to be fitted. For 
simplicity, we choose among the set of schemes satisfying (0), the one with /^s = 0, 73 
being that free parameter. For such a renormahzation scheme, 

Z^(/i) = Z(=^)(/i) exp (^^1 (S^(/^) - «^ ) ) (23) 

is, up to higher irrelevant orders, the solution to the four-loop coupled equations analogous 
to @ and (p!0D , where Z^^\jji) is the solution of the three-loop problem, with uq being 
the initial strong coupling constant at /io for both, three- and four-loops. The results of 
the fit in the energy window (2.8 -^ 4.3) GeV read (see Fig. |^) 

4^^ ^ 0.995(3) , 4M^ ^ 1.012(2) , 
Z|(a(6.8)) ^ ' ^|(«(6.8)) ^ ' 

73 = (2.2 ± 1.6) ■ 10^ [xVd.o.f. = 1.17) ; (24) 

where the errors on the ratios of Zg's do not take into account the uncertainty coming from 
the errors on the lattice spacing ratios. This uncertainty does not exceed one percent. 
When the lower limit of the energy window takes values below 2.8 GeV, the x^ /d-o.f. 
rapidly increases, which indicates the end of the four-loop matching. This is illustrated 
in Fig. ^. The upper limit is fixed by ap < 7r/2. We note also that the fitted values of 
the ratios (^^ are very close to 1.0, and somewhat larger than the results obtained by 
using the one-loop lattice perturbation theory 0. 
The estimated 73, is not as large as it might look: 

73 a/A-K , 73 a^ /{AttY , , 

^^-^ — -0.6, and — ^^^-0.05 (25) 

72 7o 

at /i = 4.0 GeV. The large error for 7J3 quoted in eq. (|^) is not surprising because the 
determination of 73 strongly depends on a^^j^(9.5 GeV). In other words, a very small 
uncertainty on oirr^ reflects in large error for 73, in our energy window (2.8 -i- 4.3) GeV. 
Thus the computation of the value for 73 from perturbation theory would considerably 
reduce the systematic uncertainty of ^^^^7^(9. 5 GeV), and hence of A'^^. 

^For comparison, the one-loop perturbative values of the two ratios are : .Z|(6.0)/2'|(6.8) = 0.929, 
and Z|(6.2)/Z|(6.8) = 0.947. in clear disagreement with the fitted values (eq. (EJ) even if the small 
uncertainty on the lattice spacing is considered. On the contrary, the fitted ratio Z|(6.2)/Z|(6.0) = 1.017 
is in good agreement with the one loop perturbative prediction, 1.019. 
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Figure 5: One-, two-, three- and four-loop perturbative curves obtained from our best 
fits of Q;j^^^(9.5GeV) and % are presented in plot (a). Notice that all these curves are 

computed with A'^*-^^ = 319 ± 14 MeV. The points correspond to the lattice evaluations 
at (3 = 6.0, 6.2 divided by our best fits of the ratios of Z^ referred to (3 = 6.8 (Z^ is taken 
to be 1 at (3 = 6.8j. Plot (b) shows the same perturbative curves and the lattice data at 
f3 = 6.8. 
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Figure 6: Plot(a): the x^/d.o.f. for the global fit obtained after the estimation of the 
four-loop contribution as a function of the lower limit for the energy window. Plot (b): 
the global fit of all the lattice results with the four-loop corrected perturbative formula 
(the full line curve on Fig. ^ from eq. ([2^. 
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At this point, we can check the consistency by performing a global fit over the entire 
window 2.8 GeV < /U < 9.5 GeV, with the values of 73 and of the ratios Z| taken from (|2^) . 
The result of such a fit is depicted in Fig. ^. The global estimate of A'^^ does not modify 
the one we obtained at large energies (5.6 GeV < /i < 9.5 GeV) given in eq. (0), whereas 
the global x^/d-o.f. is 0.79. 

Finally, we should assess the systematic uncertainty introduced by the assumption 
of the three-loop asymptoticity {i.e. 73 = 0), for 5.6 GeV < /i < 9.5 GeV. The 
simplest estimate of the non-asymptoticity errors consists in monitoring the value of A'^^, 
while varying the coefficient 73 7^ 0. I73I is reasonably bounded by Air'j^'^^ / a{n) over the 
whole window. From the previous works [0, |18| and from the present study we see that the 
non-asymptoticity has a tendency to provoke overestimates of A'^^ (the effective Aj^j^ s 
decreases as the matching is performed at higher energies). This observation is sufficient 
to exclude the negative values of 73. Then, by varying < 73 < 47r7^'^^/a(5.6GeV) ~ 
42000 we obtain the uncertainty in A^^ of ~ 20 MeV. Thus, the Landau gauge gluon 
propagator analysis results in 

^Ms _ 3^9 ^ ^4+10 ^ ^26) 

where the upper limit for the systematic uncertainty comes from the dispersion we ob- 
served by exploring the domain of three-loop good schemes, around MOM, over the (71, 72)- 
plane. It is important to note that the uncertainty on the value of A'^^ would be consider- 
ably reduced if the value 7|^o^ was known. Not only the —20 MeV would be reduced, but 
also the dispersion over the set of good schemes to four loops would be fairly restrained. 
Reciprocally, if A^^ was known accurately from any other source, we could use our 
data to fit 7^ rather accurately. For example taking A'^^ strictly equal to the central value 
in eq. (p6|), would give 



(47r)4 (47r)4 (3^ 



= 0.88 ±0.04. (27) 

AMS=319MeV 

where we have used eq. (p2D. But let us repeat, the fitted value for 7^ varies quickly when 



y^MS jg varied within the error bars, which explains the large error in (E^ or equivalently 



MOM ^MOM 

^3_^ + ^L^ = 0.9±0.6. (28) 

(47r)4 (47r)4/3o ^ ^ 

5 Discussion 

The estimate of A'^^ is significantly higher than the one obtained with the Schrodinger 
functional [|l|], A^ = 238 ± 19 MeV. Zero fiavour NRQCD results pi, although not 



directly expressed in terms of A^^, seem to agree with the result from Schrodinger func- 



tional. Estimates from string tension cover a large range of values: 244(8) MeV [|T5 
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293(18)lg3 MeV [^Oj. On the other hand, the value recently obtained directly from the 
triple gluon vertex, A^^ = 295 ±20 MeV in 0, and the less recent 340(50) MeV [§, along 
with the one obtained in this paper, favor the larger values of A'^^. The discrepancy is 
of the order of three sigma. Method based on Schrodinger functional and the one based 
on Green functions are quite different so that a direct comparison is not easy. Could it 
be that the reason for this discrepancy is simply that we did not reach a enough large 
energy ? In other words, could it be that the difference between A^^ = 238 ±19 MeV 
and the result (pB]) were simply due to the fact that next to third order terms, which are 
not used in the fit leading to (^61) , do mimic a larger A^^ ? To investigate this question 
we use a simple check: had we assumed the result for A'^^ in ref. |]18| to be the right 
asymptotic one, we would have obtained 



(47r)4 + (47r)4 p. 



= 8.42 ±0.08 , (29) 

AMS=238 MeV 



over the same energy window used for (pT]), (2.8 -^ 9.5) GeV, the x^/d-o.f. being 0.89. 
The Schrodinger functional result applied to our data would then imply that the Landau 
gauge gluon propagator is not asymptotic at three-loops at the energy scale of 9 GeV. 
The four-loop contribution in this case would be much bigger than the three-loop one 
(about four times). This seems rather unlikely! We therefore conclude that the value 
^Ms obtained by using the Schrodinger functional technique is difficult to accommodate 
with the gluon propagator data by using the three-loop expression. There may be some 
unknown systematic effect explaining this discrepancy. 

Let us finally make a comment about the convergence of the gluon propagator. The 
direct connection between the renormalization constant in the MOM-scheme and the gluon 
propagator makes the gluon momentum the natural scale in this scheme. The scales in 
the MOM-scheme are significantly larger than the corresponding ones in the MS scheme, 

typically by a factor 1/0.346 (see eq. (pll) . The Landau pole in the MOM scheme is 
around 1 GeV. As the perturbative regime of the propagator is expected to settle in when 
log(/i/A^^j^) is large enough, it is then not too surprising that high perturbative orders 
are important around ~ 3 GeV and higher energies are needed to reach a three-loop 
perturbative scaling. 

6 Conclusions 

The main goal of the present work was to go deeper into the study of the asymptoticity 
of the Landau gauge gluon propagator. The matching of its non-perturbative evaluation 
from lattice with perturbative predictions, gives us an estimate of the strong coupling 
constant and hence of A^'-^° ^ . 

We have carefully examined the lattice spacing effects, particularly the hypercubic 
artifacts, and the finite volume ones. We find a close linearity in a^pW of the gluon 
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propagator, with the slope given by eq. (|T6|) which removes efficiently the hypercubic 
artifacts. The finite volume effects, in the region of large momenta, are parametrized by 
the relation (pISl). 

After having subtracted the lattice artifacts, we found that the four-loop contribution 
is negligible above ~ 5. GeV, but becomes important below this energy, confirming the 
conclusion of ref. 0. In its turn the four-loop perturbative scaling fails below 2.8 GeV: 
the Landau gauge gluon propagator reaches very slowly the asymptotia. 

We therefore have fitted with a three-loop formula over the energy window 5.6 GeV < 
/i < 9.5 GeV. The rather good fit leads to AfJop) = 319 ± 14 MeV. A fitted four-loop 
formula has been used to extend the fit over the larger energy window (2.8 -j- 9.5) GeV. 
We have obtained a consistent description of all our lattice data. 

Our final result is 

A^ = 319 ± lAt'^ 1^^ MeV , (30) 

with the errors discussed in detail in Sec. ^. Although a combination of theoretical results 
is always delicate, we may try to combine this result with the one obtained from the study 
of the three gluon vertex |jl|, A^^ = 295 ± 20 MeV. This results in an overall fiavorless 
estimate from the gluon Green functions : 275 MeV < A^^ < 343 MeV . 
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